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Calculation of the renormalized charmed- quark mass in Lattice QCD* 

A. Bochkarev and Ph. de Forcrand 

Swiss Center for Scientific Computing, ETH-Zentrum, CH-8092 Zurich, Switzerland 

The correlator of heavy-quark currents is calculated in quenched Lattice QCD on a 16 3 x 32-lattice for /3 = 
{6,6.3}. The renormalized charmed quark mass is extracted from the short-distance part of that correlator: 
m^ fS (m c ) = f .22(5) GeV. We study the sensitivity of our data to the strong coupling constant. 



The point-to-point correlator of local heavy- 
quark currents: 

n( 9 V = * / dxe^ < 0|T{^(x)>(0)} |0 > (1) 

where = cj^c, is reliably calculable in the 
loop expansion of perturbative QCD in the vicin- 
ity of zero momentum, tf one knows the corre- 
lator from experiment or from lattice simulations 
one can then extract parameters of perturbative 
QCD such as the renormalized heavy-quark mass 
or strong coupling constant normalized on the 
heavy-quarks threshold 0. In earlier work [p|J| 
we pointed out that applicability of the pertur- 
bative loop-expansion implies that the lattice ar- 
tifacts of the short-distance part of the correla- 
tor H(q 2 ~ 0) may be studied analytically order- 
by-order in that expansion. Using this idea, our 
aim is to extract the renormalized charmed-quark 
mass m c and strong coupling constant a s (m c ). 

Following Q| we consider the ratios r n = 
M n+1 / M n of moments of the correlator (|l|): 
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to separate the short and the long distances, 
applicability of perturbative QCD near q 2 
implies the following expansion for the lower ra- 
tios r2 ) 3,4, originating from short-distances: 

r « = T^V { 1 + w »°sK) } (3) 

where {a n , uj n } are known numbers Q. The term 
oc a n comes from one loop of free charmed quarks. 
The term oc uj n comes from the two-loop diagrams 
due to one-gluon exchange. 
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Figure 1. Charmonium spectrum obtained on a 
16 3 x 32 lattice at (3 = 6 (squares) and /3 = 6.3 (dia- 
monds) and on an 8 3 x 16 lattice at /3 = 6 (crosses) 
with the clover-and-tadpole-improved action. The 
pseudoscalar mass m Vc = 2.Q7QGeV is used for nor- 
malization. The wide horizontal lines are experimen- 
tal data. 



We have calculated the correlator (|l|) in 
quenched QCD on an 8 3 x 16 lattice for f3 = 6 and 
a 16 3 x 32 lattice for (3 = {6, 6.3} with the Wilson 
tadpole-and-clover- improved fermionic action [|J , 
with 20 configurations in every case. We fixed the 
scale (lattice spacing a) from the mass of the low- 
lying resonance. To extract that mass in these rel- 
atively small volumes we incorporated the phys- 
ical continuum spectrum at high energies in the 
way it is done in the QCD sum rules approach 
to describe the behavior of the correlator at 
short distances. We used the following extension 
of the dispersion relation for the two-point corre- 
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lator II(<7 2 ) to the lattice theory 
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s + £ E M sin (<W 2 ) 



(4) 



The higher ratios rs — rii originate from long 
distances. They are saturated by the single low- 
lying resonance contribution. Our results for the 
charmonium spectrum are shown in Fig. 1. 
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Figure 2. Monte-Carlo data for the vector currents 
on a 16 3 x 32 lattice at /3 = 6 (diamonds), /3 = 15 
(squares); k = 0.1100. The lines show the ratios of 
moments for one loop of free Wilson fermions. The 
dashed line corresponds to k = 0.1100; the solid lines 
fit the 3 rd ratio of Monte-Carlo data by construction. 

We fit the lower ratios r3 4 of Monte-Carlo data 
with the free-quarks approximation to the corre- 
lator (|): 

T3 t 4 (free Wilson quarks) = r$,4 (Monte-Carlo) (5) 

The Wilson parameter k corresponding to those 
free quarks determines the rcnormalized quark 
mass. The relation between re in a hnite volume 
and the corresponding quark mass of the contin- 
uum theory was clarified in [pi. 

To specify the subtraction scheme of the renor- 
malized quark mass, one should be sensitive to 
the a s -correction in the ratios r n . If we ignore 
lattice artifacts in the coefficients w ra and take u n 
from the continuum theory pj] , choosing the MS 



subtraction scheme and a s (m c ) w 0.3 we get 
the results for the renormalized charmed-quark 
mass shown in Table 1. Since the two-loop cor- 
rection is small in the lower moments, the uncer- 
tainty in the quark mass introduced by taking u n 
from the continuum theory rather than evaluat- 
ing it on the lattice is small. 

Table 1 

The charm-quark mass obtained by fitting the ratios 
r3, r4 of the correlator on the 16 3 x 32 lattice . The 
masses m res and m c axe in units of the lattice spac- 
ing. The mass m^ is [aj/^] is in GeV with the lattice 
spacing fixed from the vector channel. Errors are sta- 
tistical. 



,3 = 6.0, re = 0.1060, (a: 


a 1.9 GeV) 


fit to T3 


Tn res 


m c 


TO*' /s [a J/v ,] 


Vc 


1.562(5) 


.618(2) 


1.208(6) 


J/l/> 


1.600(6) 


.660(1) 


1.246(6) 


Xo 


1.86(4) 


.645(2) 


1.230(8) 


fit to T4 


Tfl res 


m c 


mf s [aj/^] 


Vc 


1.562(5) 


.645(1) 


1.224(6) 


J/lp 


1.600(6) 


.677(1) 


1.258(6) 


Xo 


1.86(4) 


.686(3) 


1.27(1) 


(3 = 6.3, re = 0.1150, (at 


* 3.7 GeV) 


fit to r$ 


Tfl res 


m c 




Vc 


0.805(6) 


.307(2) 


1.149(6) 


J/lP 


0.836(7) 


.336(2) 


1.214(9) 


Xo 


0.89(6) 


.338(5) 


1.23(3) 


fit to T4 


7Tl res 


m c 


TO*' /s [a J/v ,] 


Vc 


0.805(6) 


.321(2) 


1.165(6) 


J/*l> 


0.836(7) 


.344(2) 


1.222(6) 


Xo 


0.89(6) 


.359(8) 


1.27(4) 



One can see a rather stable va lue of the 
renormalized charmed-quark mass m^ s {ra c ) = 
1.22(5)GeV, which is in good agreement with es- 
timates of the continuum theory m^f s {m c ) « 
1.23GeV. The previously reported lattice result 
is § : mp(m c ) = 1.5(3)GeV. 

We have found previously || that the lattice 
version (^) of the dispersion relation works ex- 
tremely well (i) to reproduce the lattice corre- 
lator on the one-loop level in perturbation the- 
ory, and (ii) to fit Monte-Carlo data in small 
volumes with the phenomenological form of spec- 
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Figure 3. The coefficient co n obtained via the lat- 
tice dispersion relation with the Schwinger spectral 
density on different lattices for am c = 0.6 vs. the 
continuum-theory result (dashed line). 

trum incorporating both the low-lying resonance 
and the smooth continuum spectrum at high en- 
ergies. Here we use the dispersion relation (^) 
to estimate lattice artifacts in the coefficients uj n . 
We use in eqn. (|) the Schwin ger expression for 
the spectral density, which determines the cor- 
relator of vector currents on the two-loop level. 
The quark mass is fixed in the MS'-scheme in the 
continuum theory (see jjj, where it was fixed in 
the MOM-scheme). One can see from Fig. 3 that 
the coefficient uj n obtained in this way remark- 
ably coincides with the continuum-theory results 
for the ratios of interest n — 3, 4 already on a 32 4 
lattice. Since the use of the lattice dispersion re- 
lation (W) helps to approach the continuum limit 
faster ||, Fig. 3 does suggest that one can have 
very small lattice artifacts in the coefficients u> n 
on a 32 4 -lattice for reasonable values of the quark 
mass am c ~ 0.6. Finite-size effects are seen to 
be significant on smaller lattices. 

Keeping n the same as at (3 = 6, we have gen- 
erated Monte-Carlo data on a 16 3 x 32 lattice for 
(3 = {f0, 15}. Moment ratios are shown on Fig. 2 
for = 15. The corresponding fermion mass rh c 
is very heavy, since m c • a[/3 = 15] = m c ■ a[(3 = 6]. 
Nonperturbative effects are then strongly sup- 
pressed in the range of distances probed by the 
given lattices, hence only perturbative terms sur- 



vive. Indeed one can see the Monte-Carlo data 
much better fitted with one loop of free Wilson 
fermions in a wide range of ratios at j3 = 15 
than at (3 = 6, which is in the confining phase. 
Nonetheless the deviation from the one-loop ap- 
proximation is clearly visible. Therefore, we fit 
our data with the two- loop expression (|3|), tak- 
ing the coefficients lu h from continuum theory. 
Two adjacent ratios r n , r n+ i can be used to fix 
the quark mass and the effective coupling con- 
stant a' s . Fig. 4 shows our results for a' s . Since 
af s = a s (rh c ) ~ a s (a), one can check that the rel- 
ative height of the plateaus in a' s for n = 5 — 8 
is consistent with asymptotic scaling, as expected 
at those high (3s. 
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Figure 4. Strong coupling constant for very heavy 
fermions obtained on a 16 3 x 32 lattice at = 15 
(squares) and /3 — 10 (diamonds) for n — 0.1100. 
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